Atomic multipole relaxation rates near surfaces 
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The spontaneous relaxation rates for an atom in free space and close to an absorbing surface are 
calculated to various orders of the electromagnetic multipole expansion. The spontaneous decay 
rates for dipole, quadrupole and octupole transitions are calculated in terms of their respective 
primitive electric multipole moments and the magnetic relaxation rate is calculated for the dipole 
and quadrupole transitions in terms of their respective primitive magnetic multipole moments. The 
theory of electromagnetic field quantization in magnetoelectric materials is used to derive general 
expressions for the decay rates in terms of the dyadic Green function. We focus on the decay rates 
in free space and near an infinite half space. For the decay of atoms near to an absorbing dielectric 
surface we find a hierarchy of scaling laws depending on the atom-surface distance z. 

PACS numbers: 34.35. +a, 42.50.Nn, 32.70.Jz, 32.90.+a 



I. INTRODUCTION 

The spontaneous decay of an excited atom and the 
associated release of a photon into the environment is 
due to the interaction between the atom and the vacuum 
fluctuations of the electromagnetic field. It is a funda- 
mental process that is of critical importance in the un- 
derstanding of the dynamics of free atoms. It is, there- 
fore, unsurprising that this process has been the subject 
of much study, some of which dates back almost to the 
beginnings of quantum mechanics itself [l[ . The usual ap- 
proach is to expand the atomic charge distribution into 
a sum of its multipole moments. The atom-field inter- 
action is then considered to be dominated by the linear 
coupling between the vacuum field and the dipole compo- 
nent of the charge distribution. This approximation as- 
sumes the dipole-field interaction to be sufficiently large 
that the higher-order multipole terms can be neglected. 
Such atom-dipole calculations are a mainstay of quantum 
optics textbooks (e.g. 2-5]). More recently, extensive 
work has been done on dipole interactions in more com- 
plex environments such as near absorbing surfaces @, 0] , 
near an absorbing microsphere 0, H[ and in a spherical 
microcavity Q- 

In the past few years, there has been increasing ex- 
perimental interest dipole-forbidden atomic transitions 
[3, Il0| . In particular, there has been a large body of 
work exploring the use of such higher multipole moment 
transitions in atomic clocks, because the weak nature of 
these transitions leads to a narrow line width and hence a 
much better frequency standard [ll|. Quadrupole tran- 
sitions in 199 Hg+ |H and 87 / 88 Sr+ [7jUj[| are fre _ 
quently used and even the 2 Fz — 2 Si octupole transitions 



in 171 Yb + has been studied 17]. Although there have 
been some attempts to study the theory of quadrupole 
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transitions [1814211 ] , a complete understanding of the na- 
ture of higher-order multipole decays is still lacking. 



The corresponding magnetic interaction between the 
atom and the vacuum fluctuations leads to a sponta- 
neous change of the electronic spin state. With the recent 
rise in interest in the magnetic trapping and confinement 
of atoms for quantum control experiments and quantum 
information processing [22--24J , these spin flips have be- 
come an increasingly important process. Atoms in spe- 
cific Zeeman sublevels of their hyperfine ground state can 
be magnetically trapped. However, these states are sub- 
ject to spin transitions which can take the atom from 
a trapped state to an anti-trapped state (i.e. the atom 
is actively expelled from the trap owing to the relative 
orientation of the magnetic field and atomic spin direc- 
tion). This relaxation process is a limiting factor on the 
lifetime of trapped atoms and may place serious restric- 
tions on experiments that require long trapping times 
such as those using a trapped atom as a fundamental 
qubit for quantum information processing. The dipole 
spin flip rate has been studied for a number of different 
environments, for example above conducting and super- 
conducting surfaces |25rl28| and close to metal wires [29| ■ 



In this article, we present a general theory of higher 
multipole relaxation rates. The Green function method 
for the quantization of the electromagnetic field is briefly 
discussed in Sec. [II] In Sec. IIII1 the multipole expansion 
and the multipole interaction Hamiltonian are reviewed. 
In Sec. lIVl the Heisenberg equation of motion are used to 
derive general expressions for the decay rates for various 
multipole orders from their respective interaction Hamil- 
tonians. Finally in Sees. [V] and IVI1 the decay of atoms 
in free space and close to a planar dielectric surface are 
studied in detail. Some useful but lengthy calculations, 
in particular on rotational averaging of tensors, can be 
found in the Appendices. 
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II. QUANTIZATION SCHEME 

The electromagnetic field in an absorbing magnetoelec- 
tric medium can be quantized by expanding the electric 
field in terms of the dyadic Green function (the inverse 
of the Helmholtz operator) and a set of bosonic vector 
fields, which describe the collective excitations of the field 
and the absorbing matter (for reviews, see e.g. (3fl l3ll| ) . 
The method begins with the classical Maxwell equations 
in frequency space. In the absence of free currents and 
charges these equations take the form 



V.B(r,w) = 0, 

V x E(r, uj) - iujB(r, uj) = 0, 

V-D(r,w)=0, 

V x H(r, uj) + jwD(r, uj) = 0. 



(1) 
(2) 
(3) 
(4) 



Maxwell's equations have to be supplemented by consti- 
tutive relations between the primary fields E(r,w) and 
B(r,u;), and the derived fields D(r, uj) and H(r, uj), re- 
spectively, which can be written as 



D(r,w) =e E(r,w)+P(r,w), 
H(r, uj) = K B(r, uj) - M(r, uj) 



(5) 
(6) 



[«o = 1/Mo]- The polarization field P(r, w) and the mag- 
netization field M(r,w) are, in the linear-response ap- 
proximation, related to the electric field and the magnetic 
induction by 

P(r, uj) = e [e(r, uj) - 1] E(r, uj) + P N (r, uj), (7) 
M(r,w) = /to [1 - /s(r,w)] B(r,w) - M N (r,w) (8) 

where e(r,uj) and k(t, uj) = l//i(r, uj) are the dielectric 
permittivity and the (inverse) magnetic permeability, re- 
spectively. 

Absorption in the medium is accounted for consistently 
by the addition of the noise polarization field Pfj(r,a)) 
and the noise magnetization field Mufr.w). As a result 
of this, the (frequency components of the) electric field 
obey the inhomogeneous Helmholtz equation 



V x «(r, uj) V x E(r, uj) - -^e{r, w)E(r, uj) = iwjUoj N (r, uj) 

(9) 

with the noise current density 



j N (r,w) = -JwP N (r,w) + V x M N (r,w) , 



(10) 



This equation can be formally solved using the Green 
tensor for the Helmholtz operator 

E(r,uj)=iujfio J (!VG(r,r».j N (r», (11) 

where the Green tensor G(r,r',w) has the property 

uj 2 

Vx/t(r,(j)VxG(r,r» T e(r,uj)G(r,r' ,uj) = tf(r-r') 

cr 

(12) 



Quantization of the electromagnetic field is performed 
by decomposing the noise polarization and magnetization 
fields in terms of two sets of bosonic vector fields 



P N (r,w) = £"(r,w)fe(r,w), 



M N (r,w) = iW— n"(r,uj){ m (r,uj), 



(13) 



(14) 



and imposing canonical commutation relations for them 
[A, A' = e, m] 



f A (r,a,),fl,(rV) = S xx ,S(r - t>)8(u - uj'). (15) 



Thus the frequency components of the quantized electric 
field can be written as 

E(r,«)= J d 3 r'G x (r,r',uj)-h(r',uj), (16) 

X=e 1 m 

with the magnetic induction field following from (j2Jl , and 
the abbreviations G\(r, r',uj) given by 



G e (r,r» = 
G m (r, r',uj) 



,UJ 



—e"(v',u:)G{v,v , ,uj), (17) 

C z V 7T£0 



.uj h n"{r',uj) 
C V 7T£ \fi(r',uj)\ 2 



G{t,y',uj) x ^' 



The total electric-field operator reads 



E(r) 



G?wE(r,u;) + h.c, 



(18) 



(19) 



with a similar expression holding for the induction field. 
The bosonic operators f\(r', uj) and f|(r', uj) describe col- 
lective excitations of the photonic modes and the absorb- 
ing medium and can be viewed as the generalization of 
the free-space photonic amplitude operators to arbitrary 
magnetoelectric media. The bilinear Hamiltonian 

H F = J d 3 r J dujhujil{r,uj) -ixir-u) (20) 



can be used to generate the time-dependent Maxwell 
equations from the Heisenberg equations of motion for 
the displacement field and the magnetic field. 



III. MULTIPOLAR-COUPLING HAMILTONIAN 
AND THE MULTIPOLE EXPANSION 

The macroscopic interaction between light and matter 
is commonly described using the multipolar coupling. In 
this picture, the matter is described in terms of a polar- 
ization field, which is a result of displaced charges within 
the material, and a magnetization field, which is a result 
of charge currents within the material. The interaction 
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terms take the form of linear couplings with the polar- 
ization field interacting solely with the external electric 
field and the magnetization field interacting solely with 
the external magnetic field. 

The Hamiltonian for a globally neutral system of point 
particles of mass m a and charge q a in the presence of an 
external electromagnetic field, whose centre of mass is 
at rest, has three contributions describing the medium- 
assisted quantized electromagnetic field, the free motion 
of the chargedparticles and the particle-field interaction, 
respectively (3lL l32j. 



H = Hp + H A + H int 



(21) 



with 



^=E#^ + F- / ^PaW, (22) 
^ 2m a 2e J 

a 

H int = ~Jd 3 r P A (r) ■ E(r) - J d 3 r M A (r) ■ B(r), 

(23) 

and Hp given by Eq. (|20l) . Note here that diamagnetism, 
which is quadratic in the magnetic field, is not consid- 
ered. Furthermore, since the centre of mass of the charge 
distribution is assumed to be stationary, the Rontgen 
term vanishes. By considering the divergence of the po- 
larization field and the curl of the magnetization in terms 
of displaced point charges and charge currents respec- 
tively, the polarization and magnetization fields can be 
written in integral representation as (see Appendix A) 

PaW = y2<la ds 8 [r - r A - s(r a - r A )] (r Q - r A ) , 
Jo 



(24) 



r - r A - s(r Q - r A )} 



x [(r Q - r A ) xp a -p a x (r Q - r A )} 
+ ^ 7Q S Q <5(r-f Q ). (25) 



Here S a is the spin operator for particle a located at f a , 
p a is the canonical momentum of the particle relative to 
the atomic centre of mass and 7 a is the gyromagnetic 
ratio of the particle (for an electron: j e — —eg e /(2m e ) 
with the electron g- factor g e ~ 2). Hence, the electric and 
magnetic terms in the interaction Hamiltonian become 



H 



(E) 



/ d 3 r y^gq / ds 8 [r - r A - s(f a - r A )} 

J a JO 

x (r Q - r A ) ■ E(r), (26) 



H 



(B) 



r - r A - s(r a - r A )] 



x [(r Q - r A ) xp„-p a x (r Q - r A )] • B(r) 

- f d 3 r J2 J a S a 6 (t ~ fa) ' B(r). (27) 

A. Electric multipole expansion 

The interaction Hamiltonians and (|2~T1) arc diffi- 
cult to handle and some simplifications are necessary to 
proceed. The usual approach is to expand the polariza- 
tion and magnetization fields in terms of their multipole 
moments. The procedure is outlined here. More infor- 
mation on the both the classical and quantum multipole 
expansion can be found in Refs. [33j-|35j. 

Consider first the electric term of the interaction 
Hamiltonian (|26[) . The 8- function can be expanded about 
the point r = r A 

^ff = "£<?«/ d\ j ds (r a - v A ) ■ E(r) 
J Jo 

x \ 8(r- t a ) - s(r a - r A ) ■ V8{r - r A ) 



^[(f Q - v A ) ■ V} 2 S(r - v A ) + 0[(r a - r A f] 



(28) 



Both integrals can now be performed. The spatial inte- 
gral is evaluated by integration by parts on the deriva- 
tives of the 8- function. As the vector (r a — r A ) is not a 
function of the spatial variable r, the derivatives act only 
on the electric field. Thus the expansion becomes 



H, 



(E) 



q a (r a - r A ) ■ ~E(r A ) 



q a [(r a - r A ) ® (f Q - r A )] • [V <3 E(r) 



2! ^ 

a 

~ 3? £ qa ^ a - t a) ® (i a - ?a) ® (f a - r A )] 

a 

• [V ® V ® E(r)]| r=rA + O [(r Q - r A ) 4 ] , (29) 

where the symbol • denotes the Hadamard product (A • 
B = Ai 1 ...i n Bi 1 ...i n ). Equation (|2"9")l can be written in 
a simpler form by defining multipole moment (tensor) 
operators. The primitive electric dipole, quadrupole and 
octupolc moment operators are defined as 



d = 2_^q a (r a - ta), 



(30) 



d (4) = ^£g Q (f Q -r A )«,(f a -r j4 ), (31) 

a 

d (8) = -y q a (r a - r A ) ® (r Q - r A ) ® (r a - r A ), 



(32) 
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respectively. Note here that the multipole moments are 
defined as in Refs. 0,[35[ with the coefficients of the Tay- 
lor expansion absorbed into the definition of the moment. 
As a result, the electric multipole interaction Hamilto- 
nian can be written as 



H, 



(E) 



d-E{r A ) ~d (4 ^ *[V > 
d (8) • [V ® V ® E(r)] 



E(r)] | r=rA 

r=r A +0 [(T a 



(33) 



B. Magnetic multipole expansion 



A similar expansion can be performed for the magnetic 
term of the interaction Hamiltonian (|27|) . Once again the 
(^-function can be expanded about the point r = r A ■ Per- 
forming the integrals as before and defining the primitive 
magnetic dipole and quadrupole moments as 



(34) 



a 

3! ^-^ m a I 

a v 

+ ((r a - r A ) <g> [(r Q - r A ) x p a ] ) T | 

+ ^2(r a - r A ) (X>7 a S Q , (35) 



respectively, leads to a magnetic multipole interaction 
Hamiltonian of the form 



H$ = - m • B(r A ) - mW . [V ® B(r)]| r= 
- O [(r a - r A ) 3 ] . 



(36) 



IV. ATOMIC EQUATIONS OF MOTION 

The dynamics of atoms subject to electric dipole in- 
teraction with the medium-assisted electromagnetic field 
has been studied extensively in the literature p-Q ■ More 
recently, atoms subject to electric quadrupole interac- 
tions [2l[ and magnetic dipole interactions [26l - l29j have 
been considered. In order to analyse the effect of higher- 
order multipole interactions on the atomic decay rates, 
we derive the equation of motion for an atom subject to a 
single, specific higher-order interaction. We develop the 
Heisenberg equation of motion for general multipoles and 
derive a general expression for multipole driven relax- 
ation rates. It is assumed that each transition is driven 
by a single multipole interaction term and hence each 
multipole order can be studied independently. 

Beginning with Eq. (|21l) . the microscopic multipolar- 
coupling Hamiltonian can be written as 



H = Hp 



Hi 



int • 



(37) 



Here an — \i)(i\ are the projection operators onto the 
energy eigenstates \i) with energies frioi of the atomic 



Hamiltonian, and Hp is given in Eq. (|20|) . We define 
an operator Q r that corresponds to the (differential) op- 
erator acting on the electric field in the relevant term 
of the multipole expansion [e.g. for the electric dipole 
Q r = (d-), for the electric quadrupole Q r = (cT 4 ) • V r <8>), 
etc]. Furthermore, we can rewrite the interaction Hamil- 
tonian as 

flint - Q r E(r) | r=rA 



= -E 2 ^E(r) \r=r A , 



*jj I E ( r ) i*=*a 



(38) 



where the electric field is given by Eqs. ([Toll and (fH?|) . 
and the completeness relation an — 1 has been used. 
The Heisenberg equations of motion for the atomic flip 
operators and the bosonic field operators are thus 

f A (r,w) = -iu)f x (r,u>) - i^gj. y (r',r,w)<7.y (39) 



and 



j&ij -i J2 d3s 



duj 



A — e , in k 

[gxjk (r, s, uj)a ik - gx,ki (r, s, u)a k j] ■ f\(s, w) + h.c. 

(40) 

where the coupling tensors gx ij(T,r' ,u) are defined by 



gx,ij{r,r',u)) = — Q£,G A (r, r', w 
in 



(41) 



Our goal is to study the dynamics of the atom under 
the influence of an external field, hence the next step is 
to remove the electromagnetic degrees of freedom. This 
is done by solving the equation of motion for the field 
(131)1) and substituting the solution into the equation of 
motion of the atom (|4"0")l . Equation (|3"9"|) can be solved in 
the Markov approximation which assumes that atom to 
be memoryless and hence not affected by previous field 
states. Writing the atomic operator as a slowly varying 
envelope function and a term oscillating with the atomic 
transition frequency, a^ — cr^e*"'"*, Eq. (|39|) can be in- 
tegrated to give 

f A (r jW ) = e— *fI reo (r, W )+^ g ^(r',r, W )C( W -^)a lJ 

y 

where £(w — u>ij), in the long time limit, is defined as 



C(w - Uij) = V 



+ iTrS(oj — ujij) (43) 



with V denoting the principal value. Resubstituting the 
expression for the field into the equation of motion for 
the atom gives 
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% = iWij&ij - i ^2 ^2jd 3 s J due lut [gxjk(r, s, w)<x ii: - a\,fci(r, s, w)<j fci ] • f| roc (s, w) 

A— e,m fe 

X! E d 3 s dLd[g x ^k(r,s,u)a ik - g\M(r,s,u)a k j] ■ gl m (r' ,s,u)((u} - u} pq )a pq . (44) 



A— e,m fc,p,g * 



The first term in Eq. (j44|) corresponds to the free mo- 
tion of the atom, the second to the interaction of the 
atom with the free external field. The third term cor- 
responds to an atom self interaction mediated by the 
medium- assisted electromagnetic field. It is this term 
that leads to spontaneous relaxation processes. 

State populations are described by the projection op- 
erator a a. Since the free motion of a projection operator 
on to a state is static the first term in Eq. (|44| van- 
ishes. If there is no external field then the second term 
also vanishes. The problem can be simplified further by 
studying only the two relevant atomic levels that are in- 
volved in the transition. Since the multipole moment 
matrix element contained in Q£- vanishes for paired in- 
dices, gu(r,s,ui) vanishes and the evolution of the excited 
state (here labelled as 2) can be reduced to 



<5"22 = i ^2 j ' J dtu 



g A42 (r,s,cj) ■gl 21 (r',s,u))a n ((u - w 2 i) 
■ 9A,2i(r,s, w) g* x 12 (r', s, w)<7 22 C(u; - w i2 



The state labelled 1 indicates the ground state. Solving 
this equation for the excited-state population gives 



022 (*) 



_ -i(a-ir)t 



0-22(0) +i ^2 J dt ' J d3s J duj 



X—e.m 



x g\,i2(r,s,u) ■ gl 21 {r' ,s,uj)an(t')C(uj - uj 2 i) 



(46) 



Here Cu = u + 5uj is the shifted atomic transition fre- 
quency [3l|. In the absence of an external field the level 
shift is small and hence can be neglected. The relaxation 
rate appears as an exponential decay factor in the free 
evolution of the excited-state projection operator and is 
given by 

r = 7r / d 3 s9A,2i(r,s,w 2 i) ■ g* XA2 (r',s,u 2 i). (47) 

\ — « 



Expanding the coupling tensors gx,ij (r, s, uS) and using 



the integral relation for the Green functions, 
/d 3 sG A (r,s,w)-G A (s,r» 

\ — ** 



hio 2 

7TC 2 £o 



ImG(r,r», (48) 



the final expression for the decay rate is found to be 



J 21 



Ql x lmG{v,v',u>2i)Ql 2 



(49) 



Equation (l4"!)f is a general result and hence true for all 
multipole orders. In order to obtain any relaxation rate 
one only needs to determine the form of the differential 
operator Q£- for the relevant order, which itself can be 
easily found from the multipole expansions (|33p and (|36[) 
of the electromagnetic field. Together with the definition 
of the multipole differential operator Q^- , Eq. (|49]) con- 
stitutes the main result of this article. One should note 
that Eq. (|4"9")l is consistant with the results of Ref. [2l[ 
where a general expression for the relaxation rate was 
derived using a classical approach. Substituting in the 
relevant expressions for Q\, the decay rates for various 



(45) multipole order can be found: 



r = d 2 i • lmG{r A ,r A ,U2i) ■ d 12 , 

,2 



he 2 Eg 

lim ^-d ( 2 ? • V®ImG(r,r',u;2i) 

r,r'-s-r,4 hc z £() 



(50) 
(51) 



T = lim — — 



d$ • ( V <8> V) ® Im G(r, r', w 2 i) <8 (V <8 • d[ 



(8) 
2 ! 



r = lim ^m 2 i • V x Im G(r, r', u 2 i) x V' ■ m 12 , 



(52) 
(53) 



r= lim 

r,r'— *ta h 



V X ImG(r,r',w 2 i) X V' Ov'-rn^ 



.(*) 



(54) 



In these expressions, Eqs. (1501) — (|52|) are the electric 
dipole, quadrupole, and octupole transition rates, respec- 
tively, and Eqs. (|53l) and (|54p the magnetic dipole and 
quadrupole transition rates. 
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V. MULTIPOLE DECAY RATES IN FREE 
SPACE 

The expressions for the spontaneous decay and mag- 
netic spin flip rates in Eqs. ([50)1 — (|54p are valid for any 
environment with the geometry of the system contained 
in the, as yet unspecified, Green function. Calculating 
rates for different environments requires knowledge of the 
specific Green function for that system. For complicated 
geometries its calculation can be highly involved, but in 
some simple cases it is analytically known. In this sec- 
tion, the Green function for free space is used to calculate 
the vacuum decay rates. In the next section, the Green 
function for an infinite half space will be used to calcu- 
late the decay rates of atoms near an absorbing dielectric 
surface. 

The Green function for free space, in component form, 
is given by 

Gg(r,r» = (ffig + g^) J^'-^l ' (55) 

where q = lo/c is the wave number and the greek indices 
run over the Cartesian coordinates x, y and z. 



free space Green function for coincident spatial variables 
are easily found. Note that the both the Green func- 
tion and are symmetric in their respective indices. 
Using these symmetries it is possible to show that the 
only non- vanishing terms are those with paired indices in 
the outer product of the two quadrupole moment vectors 
(sometimes referred to as the quadrupole polarizability 
tensor). After collecting terms, the quadrupole sponta- 
neous decay rate in free space is found to be 

Note here that the primitive moments have been used. 
In certain circumstances (e.g. when the external field is 
divergence free), the quadrupole moment tensor can be 
made traceless [33|, [3J| . In these cases the second term 
in Eq. (jST))) vanishes. This returns the result which is 
consistent with Ref. [2l| where the traceless quadrupole 
moment tensor is used. However, traceless moments are 
not valid for all physical situations (for example, if there 
are other sources in the external field) and hence must be 
used with care. In contrast, the primitive moments are 
always valid and hence the full form of Eq. (jBTJ)) is true 
for all physical situations. 



1. Electric dipole decay rate 

Extensive work has been done on the dipole transition 
For completeness the result is stated here. The 
spontaneous decay rate for an electric dipole transition 
is given in component form by 



r = lim 



d a ^ImG a ^(r,r». (56) 



Here the state indices have been dropped for clarity 
(hence o->i2 = u>2i — w and di2 = d2i = d). From 
Eq. (|55p. the imaginary part of the free space Green func- 
tion in the coincidence limit r — > r' is found to be 



lm G { ° } p (r A ,r A ,uj) = —5 af3 . 



(57) 



Hence the free space dipole spontaneous decay rate is 
given by the well known formula 



To 



67rhc 3 So 



dad a 



(58) 



2. Electric quadrupole decay rate 

The quadrupole spontaneous decay rate in component 
form is given by 

2 

r= lim -^dWd$d a d'ImGp S (r,r',u). (59) 

Once again the state indices have been dropped for clar- 
ity. The double derivatives of the imaginary part of the 



3. Electric octupole decay rate 

The octupole spontaneous decay rate in component 
form is given by 



U) 



The derivatives of the imaginary part of this function for 
coincident spatial variables are again easily found. Here 
d^l is symmetric in all its indices. Thus, again, all 
non- vanishing terms have paired indices in the octupole 
polarizability tensor. Hence the octupole spontaneous 
decay rate in free space can be written as 



Tn = 



2w 7 



I057T?lC 7 £o 



(62) 



As for the quadrupole decay rate, the primitive octupole 
moments have been used. In certain circumstances trace 
constraints cause the second term in Eq. (|62p to vanish. 
Note that for successive multipole orders the free space 
decay rates change by a factor proportional to lo 2 R 2 /c 2 
where R is the characteristic length of the charge distri- 
bution of the atom. 



4- Magnetic dipole decay rate 

Let us now turn to magnetic transitions. The calcu- 
lations are slightly more involved due to the presence of 
the spatial derivatives acting upon the Green tensor. In 
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component form the magnetic dipole transition rate is 
given by 

r = lim ^m a m l j,e a j3 7 e[ lv \d0d' J/ lmGy\(r,r',Lj). 

r,r'—tTA h 

(63) 

Calculating the derivatives one finds that the double curl 
of the Green function in the coincidence limit is 

3 

lim £^ 1 £ F ja^>G;, ) (r,r» = - — =6 a p. (64) 

r,r'->r A ' D7TC 

Hence the magnetic dipole spin flip rate is given by 
1 o = „ z ^ m a m c 



6Trhc 3 



(65) 



This is consistent with previous calculations based on 
dyadic Green functions [29| as well as on employing 
Fermi's Golden Rule El. 



5. Magnetic quadrupole decay rate 

Finally, we calculate the quadrupole moment spin flip 
rate. The notable complication at this order, compared 
to that of the electric quadrupole moment, is that the 
primitive magnetic quadrupole moment is not symmetric 
in its indices. In component form the quadrupole spin flip 
rate is given by 

r,r'— >r A fo * 

x3 o a 7 ^9iImG fc (r,r». (66) 

Computing the derivatives one finds the magnetic 
quadrupole decay rate in free space to be 



(4) (4) 1 / (4) (4) . ( 4 ) (4)\\ 



15irhc 5 

(67) 

As with the electric decay rates, successive magnetic mul- 
tipole orders change by a factor proportional to lj 2 R 2 /c 2 . 
From the above recipe, multipole decay rates of even 
higher order can be easily derived, if needed. 



VI. MULTIPOLE DECAY RATES NEAR 
ABSORBING SURFACES 

In this section the decay rates of atoms near an absorb- 
ing surface are calculated. We shall consider surfaces that 
are not bianisotropic (i.e. no cross-responses between 
electric and magnetic degrees of freedom) . Note also that 
diamagnetic materials are excluded from this considera- 
tion as diamagnetism is an inherently nonlinear process. 
Furthermore, we are interested in atoms that are sub- 
ject to either an electric or magnetic interaction. Thus 
there are in principle four possible combinations of elec- 
tric/magnetic atoms interacting with electric/magnetic 



surfaces. However, the work involved can be reduced 
when one considers the duality properties of the elec- 
tric and magnetic fields. It was Silberstein [13] who first 
suggested the existence of symmetries in Maxwell's equa- 
tions under exchange of the electric and ma gne tic fields. 
This idea was further developed in Refs. [38l.l39j. and the 
conce pt h as been more formally established recently in 
Refs. [40|, |4l| in the context of macroscopic QED. It has 
been shown that under global exchange of the two fields 
various QED quantities, in particular relaxation rates, 
remain unchanged. 

Hence, in order to calculate the magnetic relaxation 
rate for an atom above a purely magnetic surface one 
would only need to calculate the electric relaxation rate 
for an atom near a purely electric surface and apply the 
global transformations e — > /i and d — > — rh/c. Similarly, 
in order to find the electric relaxation rate of an atom 
above a purely magnetic surface one only needs to cal- 
culate the magnetic relaxation rate for an atom near a 
purely electric surface and then apply the same global 
transformation. In light of this duality, we will only cal- 
culate two of the four combinations and then deduce the 
other two from duality. We shall formally calculate the 
interaction of atoms above a purely electric surface since 
it is these systems, owing to the dominant strength of 
electric interactions, that are of most practical interest. 

In the following, we envisage a medium that is in- 
finitely extended along the (x, y)-directions and layered 
in the z direction such that 



e(r,u>) 



1, z > 
e(uj), z < 



(68) 



that is, an absorbing dielectric material of permittivity 
s(u>) occupying the lower half space with a single interface 
to free space at z = 0. The Green function can be split 
into three parts, 

G(r,r» 

G bulk (r,r',oj) +R(r,r',u), z,z' > or z, z' < 
T(r, r',uj), otherwise 

(69) 

where T(r,r',u;) is the Green function for waves 
transmitted through the interface, U(r, r',uj) is the 
Green function for waves reflected at the interface and 
G bulk (r, r',uj) is the Green function for a infinitely homo- 
geneous, isotropic (bulk) medium of permittivity equal to 
that of the layer where the spatial variables are located. 
From the general expression for the decay rate (|4"9"| it is 
clear that the Green function has to be evaluated in the 
limit of coinciding spatial arguments. Hence, the trans- 
mission part T(r,r',cj) does not play any role. Further- 
more, since the atom is located in free space above the 
dielectric surface, the bulk contribution reduces to the 
free-space Green function (f5"5"|) . Hence the decay rate of 
an atom close to an absorbing surface can be decomposed 
into the free space decay rate and a reflective correction 
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owing to the presence of the surface 

r = r + r (s \ 



(70) 



The reflective part of the Green function can be written 
as the partial Fourier transform of its components 



R{t,t',u) 



d 2 k 



(2nf 



li(k||,z,z'a;)e 2k l|- (r 



(71) 



where rii = (x,y,0) and k|| = (k x ,k y ,0) are vectors re- 
stricted to the (x, y)-plane parallel to the dielectric inter- 
face, and A; 1 1 — |k|||. The Green tensor components for 
a dielectric material in this specific geometry are listed 
in Appendix B. The symbols r™ and r TE denote the 
Fresnel coefficients for TM and TE waves 

e(uj)(3 z> o - (3 Z<0 jte _ Pz>o — Pz<o 



TM 



e(ui)P z>0 + P z< o' 



Pz>0 + Pz<0 



with 



Pz<o = \Jq 2 e{u) - fcjj, f3 z>0 = ^j~q 2 



(73) 



where Im/3 2 > and q = lu/c. The generic form of 
the Green function does not lend itself to analytical in- 
vestigations. However, in certain regimes it is possible to 
simplify the expression greatly. Here we consider the two 
limiting cases of near-field and far-field regimes. 



A. Near field 

For a purely dielectric material the near-field ap- 
proximation assumes that the atom-surface distance 
is less than the effective transition wavelength, thus 
z\y/e(u)\u/c <C 1 and hence q 2 , q 2 \s(uj)\ <C ftjj. Thus it 
is possible to expand the expressions for (3 z <o and Pz>o 
about q 2 /k 2 e(u>) and q 2 /k 2 , respectively. As a result, 



z<0 



ik\ 



P. 



z>0 



ik 




i 



(74) 



By considering only the leading order terms in these ex- 
pansions, we can find the leading order correction to the 
decay rates associated with electric transitions. Thus, 



11 £{UJ) + 1 



„TE 



(75) 

and hence the reflective part of the Green function, in 
the near-field approximation, to leading order is 



-R(kii , z, z'lu) 



1 



e(w) + 1 



T7 fc„ ~ lK x 



-fc||(W + l 2 'l) 



\ ik x 



k l{ iky 



(76) 



For magnetic transitions near purely dielectric bodies, 
to leading order in /3 2 <o and /3z>o there is no correction 
to the decay rate. In this case, it is necessary to include 
the next-to-leading order terms in /3 2 <o and P z >o, 



j3 z<0 = ik\\-i—e(u), f3 z>0 



ik\\ 



2kw 



Evaluating the double curl [V x R(kii , z, z'uj) x V] , ex- 
panding all occurrences of /3 2 <o and /3 2 >o as above and 
neglecting terms of order O [q A /k 2 ^j one finds that terms 
proportional to r™ vanish. Only terms proportional to 



„TE 



remain which itself can be approximated as 



TE 



p. 



z>0 



P 



z<0 



P. 



z>0 



P 



z<0 



Hence, the double curl of the reflective part of the Green 
function, in the near field approximation, to next-to- 
leading order is 



[V x R(k h z } z'uj) x V] 



[e(u) 1] 



,-k|i(H+|z'|) 



k x ky 



k x ky 



%k x k 1 1 



% k x k || % ky k 1 1 



(79) 



Equation (f76| will be used, together with Eqs. ([50)) 
([52]). to find the transition rates associated with electric 
transitions of an atom above a purely dielectric surface 
and Eq. ([79]) will be used, together with Eqs. (|53|) and 
(1541 . to calculate the corresponding transition rates as- 
sociated with magnetic transitions. 



1. Electric multipole transitions 

In the following, we specify our general near-field re- 
sults to the electric dipole, quadrupole and octupole tran- 
sitions. We start by recalling the well-known results for 
electric-dipole transitions. Upon using Eq. (f56|). the re- 
flective correction to the spontaneous dipole decay rate 
is found to be \£MM 



r (S) 



1 



16irhe a z 3 \e(cj) + 1| 2 



{d 2 x + d 2 y + 2d 2 z }. (80) 



Although this is a general result, in many cases such de- 
tail is not necessary. The above result assumes knowl- 
edge of the individual components of the dipole moment 
vector, and hence the orientation of the dipole moment. 

For a free atom ensemble the dipole moment orienta- 
tion of each individual atom is generally unknown. For 
sufficiently large ensembles of atoms, or for cases where 
the dipole orientation is unknown, it is possible to rota- 
tionally average the tensor formed by the outer product 
of the two dipole moment vectors (also known as the 
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polarizability tensor) over the sphere to obtain an ex- 
pression for the decay rate that is isotropic with respect 
to the dipole moment orientation. The method used for 
the rotational averaging involves contracting the polar- 
izability tensor with an averaging tensor (see Appendix 
C), 



I?fd a dp. 



3 'V 

Hence the spontaneous dipole decay rate becomes 
,2 



hc 2 e ' I if dad l 3lTI1 Gij{*A,*A,U). 



The rank-2 averaging tensor is 



(81) 



(82) 



rf (4) .(4) 




j(4) ,(4) 


A n e 








3 


d (4, d (4) 
'-ijyy *-*>yy 


3 


d (4, d (4) 


8 


d (4) d (4) 

^•xy u>xy 


1 


d (4) d (4) 


4 


d (4) d (4) 


4 


d (4) d (4) 

i^xy , -*>yx 


1 


d (4) d (4) 

U/y z ^zy 


4 


d (4) d (4) 


4 


d (4) d (4) 

i^yx t^xy 


1 


d (4) d (4) 

U>g,y Uiyz 


4 


d (4) d (4) 


4 


d (4) d (4) 


1 


d (4) d (4) 


4 


d (4) d (4) 


4 


d (4) d (4) 

l^xx ^yy 


1 


d (4) d (4) 


-4 


d (4) d (4) 


-4 


d (4) d (4) 

yy ^xx 


1 


d (4) d (4) 
u zz ^yy 


-4 




-4 



TABLE I: Table of coefficients for the quadrupole spontaneous 
decay rate of an atom close to an absorbing surface. 



which leads to a decay rate of 
,2 



r. — 

x iso — 



0J- 



3hc 2 EQ 



d a d a Im Gu(r A ,r A ,Lj). 



(83) 



(84) 



Note that the free-space decay rate is unchanged since the 
free-space Green function for coincident spatial variables 
is already spherically symmetric. The reflective correc- 
tion to the decay rate thus becomes 



r 



1 



6 » 



127rfe z 3 |e(w) + 1| 

e"(w) / r N 3 



2 dfyido. 



( c 



2|e(w) + 1| 2 Vwz 



r , 



(85) 



where To is the dipole spontaneous decay rate (|58l) in free 
space. 

Similar considerations lead to the decay rates associ- 
ated with higher-order electric multipole transitions. The 
reflective correction to the quadrupole spontaneous decay 
is given by [recall Eq. ([591 ] 



r (S) 



—d [i) d {i) 



Una 



x dad'^lrn 



d 2 k\\ 



ik l|-( r H- r N 



(86) 



Remembering that the quadrupole moment tensor is 
symmetric in its indices gives a reflective correction to 
the quadrupole decay rate of 



r(s) 



64irhs z 5 \e(uj) + 



W {A a ^/;]}, (87) 



where the term in curly brackets is a sum over the 
quadrupole moment tensor components. The coefficients 
A a j3~f8 for the non- vanishing tensor components are given 
in Table HI We can again rotationally average over the 
quadrupole polarizability tensor to obtain an isotropic 



expression for the decay rate of an unorientated ensem- 
ble 



r iso = lim 
with 



,(4)iso ,(4)iso _ ja^S ,(4) ,(4) 
u ij a km ~ 1 ijkm u al3 u 'jS ■ 



(89) 



The full expression for 1°^.^ is given in Appendix C. We 
thus find a reflective correction to the spherically aver- 
aged quadrupole decay rate of 



3 g "M L {4) d (4) - -d^d [i) 
lOTrfeo^ 5 Hu>) + 1| 2 \ d <*e d <*P 3 aa ^ 

s"{uj) 



R g W (JLYr 

|eH + l| 2 \uz) °' 



(90) 



where To is the quadrupole spontaneous decay rate (fST)]) 
in free space. 

The same procedure can be repeated for the electric 
octupole transition whose decay rate is given by [cf. 

Eq. HEED] 



x Im 



d 2 k\\ 



(91) 



By the same token, the reflective correction to the oc- 
tupole decay rate is 



r (s) = 



45 



2567rfe z 7 |e(w) + 1| 



f A j(8) ,(8) \ 



(92) 

with coefficients Aap^s/iv given in Table [TT1 Note that the 
coefficients are unchanged under all index permutations 
of the octupole moment tensor and by the commutation 
of the two tensors. For unorientated ensembles the decay 
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d l8) .(8) 



d {8) d {i) 

U 'XXX U >XXX 

U-yyyU-yyy 
d (8) d (8) 

^ZZZ^ZZZ 

dxxydxxy ~\~ P&V7YI. 
dxxzdxxz ~\~ PGTTTI. 

dyyxdyyx ~\~ P&TTTl. 
dyyzdyyz H~~ PGV7TI. 

dfLdfL + perm. 
di%dzJy +perm. 

dxyydxxx ~\~ P&V77X. 



-jm — 



5 
5 

16 
1 
6 
1 
6 
8 
8 
1 



d (8) d (8) 



d (8) d (8) 

Ui yxx u >yyy 

d (8) d (8) 

u >yzz u 'yyy 

d (8) d {8) 

Uj zxx LL zzz 

d (8) d (8) 

u, zyy u, zzz 

d (8) d (8) 

u 'xyy Ui xzz 

d {8) d {8) 

^yzz^yxx 

d (8) d (8) 

u 'Zxx Uj zyy 

d (8) d (8) 

^xyz^xyz 



+ perm. 
+ perm. 
+ perm. 
+ perm. 
+ perm. 
+ perm. 
+ perm. 
+ perm. 
+ perm. 



715) 



-6 
1 

-6 
-8 
-8 
-2 
-2 
2 
2 



TABLE II: Table of coefficients for the octupole spontaneous 
decay rate of an atom close to an absorbing surface. 



rate becomes 

r = lim 
with 



v- „ '—dm °<# s °di W;im G fcr (r, r», 

r,r'->rj ftC £q j 

(93) 



j(8)»so ,(8)»so _ ja^SriCm ,(8) 
"ijfc "pqr _ 1 ijkpqr u al3f U 8riC' 



(94) 



The full expression for 1°^,^ is given in Appendix C. 
The reflective correction thus becomes 

(S) _ 9 g "H f d (8) ,(8) _ 1 ,(8) ,(8) 



135 e"(uj) 
~2~|e(u;) + l| 2 W 



— r 



0. 



(95) 



where To is the octupole spontaneous decay rate (|62p in 
free space. 

Equations ([55]) , (ji?0"[) and (|95l) reveal that with increas- 
ing order of the electric multipole, the near-field scaling 
law of the relevant spontaneous decay rates change by a 
factor proportional to R 2 /z 2 where R is the character- 
istic charge separation length of the atom and z is the 
atom-surface distance. 



2. Magnetic multipole transitions 

We complete this study of near-field decay rates by 
investigating magnetic multipole transitions. Magnetic 
dipole transitions in the presence of dielectric materi- 
als have been intensively studied previously pBI . [2fj| , and 
their results are presented here for completeness. The re- 
flective correction to the magnetic dipole transition rate 



is given by [recall Eq. ([63 

r (S) 

r,r'- 

d 2 k 



lim —m a mp 

r,r'— >ta h 



x Im 



/4[^%,y',)xV] / %<'- 



4) 



(96) 



which, after inserting the expressions for the Green ten- 
sor, leads to 



-i(S) _ Mo^~ 



;£ »{r 



2m 2 }. 



(97) 



This is consistent with parts of the results in Refs. [25l . 
|26| . corresponding to spin flip rates for an atom above 
a metallic film in the case of large skin depth and film 
thickness. 

In the case of electric multipole transitions, it made 
sense to consider unoriented atomic ensembles and to 
define rotationally averaged spontaneous decay rates. In 
principle the same averaging process could be performed 
here. However, the result would not be physically mean- 
ingful since the atomic spin is quantized about a specific 
quantization axis and hence cannot have an arbitrary di- 
rection. Furthermore, in most practical situations the 
spins of an ensemble of atoms are aligned by external 
magnetic fields and hence assume a particular orienta- 
tion. Thus the rotationally averaged quantity is only of 
academic interest and hence its calculation is renounced. 

Finally, the reflective contribution to the magnetic 
quadrupole decay rate is given by [recall Eq. (|6"6"])] 



d 2 k 



x Im 



/^[ VXjR(k "' Z ' Z ' w)xV ]^ elk "' (r "" rI ' ; 



which, in the near- field limit, can be written as 



r (5) = 



£"{u){A aMS w^lm^}} 



5l2Trhc 2 z 3 



(98) 



(99) 



The coefficients for the non- vanishing components of the 
tensor A a /3-yS are given in Table Q] with the obvious in- 
terchange dij O rriij. Note here that the coefficients are 
the same as those for the electric quadrupole decay rate. 
This is because they are properties of the geometry of the 
system and hence come from the Green function. For the 
same reason as given for the dipole case rotational aver- 
aging is not performed. As with the spontaneous decay 
rates, the magnetic transition rates also change by a fac- 
tor proportional to R 2 /z 2 . 



B. Far field 

In the far-field limit one assumes that the distance to 
the surface z is large compared to the effective transition 
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wavelength and thus z\yj e(u>)\ui / c ;§> 1. In that limit one 
can apply the method of stationary phase to compute the 
first-order contribution to the Fourier integral ([71]) . The 
integral can be split into contributions from propagating 
waves (fcii < q) and evanescent waves (jfeii > q), the latter 
of which can be neglected in the far field. The propagat- 
ing part consists of a product of an amplitude function 
with a sinusoidal waveform whose frequency is modu- 
lated by a further function with one stationary point in 
the range of integration. At the stationary point fc|| = 
we have 



Pz<o = qV £ ( UJ ) > Pz>o = q, 



and therefore 



r TE ^ 1 \/ £ ( UJ ) r TM ^ j,TE 

l + v^H' 



(100) 



(101) 



Thus, in the far-field limit, the reflective part of the Green 
function becomes 



R 



n I \ i TE i \/l 2 - k «(\ z \ + \ z '\) j- ft -i n\ 

(km , z, z u>) = — r e V « diag(l,l,0) 



2q 



(102) 

Analogously, the double curl of the reflective part of the 
Green function in the far-field limit is found to be 

[V x R(k||,z,z'cj) x V] 

= !J r ^ e 4 V^f (|2|+|2,|) diag(l, 1, 0) . (103) 



respectively. 

As with the near-field decay rates these results can be 
rotationally averaged using the appropriate rotational av- 
eraging tensors. For an unorientated ensemble of atoms 
the Ith electric multipole decay rate becomes 



r 



(S} -(-D l+1 ^\r TE \ 



2ujz 



argr 



TE 



(-) 



(107) 

where Tq is the respective spontaneous decay rate in free 
space ([55]). (|60)) or (|62j) for dipole, quadrupole or octupole 
transitions, respectively. We see that successive multi- 
pole orders change by a factor proportional to lj 2 R 2 /c 2 , 
just as the free space ones do. 



2. Magnetic multipole transition rates 

The reflective corrections to the transition rates asso- 
ciated with magnetic dipole or quadrupole transitions is 
given by Eqs. (1631) and (j66]), respectively. By perform- 
ing the inverse Fourier transform and taking the limit as 
r, r > ta, the reflective correction to the dipole spin flip 
rate is found to be 

T-i(S) MO^ 2 TE I • f^LOZ Te\ r 2 . 2l 

l w =— - — r I sin I hargr \{m x +m \, 



87T He 2 2 



and for the quadrupole rate one finds analogously 



v J 
(108) 



1. Electric multipole transition rates 

The reflective corrections to the far-field electric dipole, 
quadrupole and octupole transition rates are given by 
Eqs. (|56"|) . (15^1) and (|6"Tj) . respectively. By performing the 
inverse Fourier transform over feii for each component and 
taking the limit as r,r' — > ta, the reflective correction 
to the spontaneous decay rates associated with dipole 
transitions is found to be 

,2 



87rSc 2 £Q' 



\r TE \sm 



IbJZ 



argr^j {dj+dj}. 

(104) 

Note here that r TE is the approximated coefficient given 
by Eq. (|101[) . Similar considerations lead to the far-field 
contributions to the decay rate associated with electric 
quadrupole transitions, 



<s) _ 



87T?IC 4 £0' 



I raisin 



2ujz 



argr 



TE 



{4t ) 4t ) + < ) 4t ) } (io5) 



and octupole transitions, 
, .6 

r (S) = 



8Trhc e eoZ 



r TB |sin 



— + argr /A 



X ^zzx^zzx + 4zy4zj/| ' (106) 



r (s) = 



\r TE \sm. 



2ujz 
c 



8irhc 4 z 
x{™i 4 Jmi 4 i+m( 4 >W}. 



argr 



TE 



(109) 



As with the near-field transition rates the presence of 
a fixed quantization axis prevents rotational averag- 
ing. Once again we see that successive multipole orders 
change by a factor proportional to ui 2 R 2 /c 2 . 

In the far-field regime, both the electric and magnetic 
multipole rates show a 1/z-behaviour. This is to be ex- 
pected as the far fields created by all multipoles have the 
same distance scaling. 



C. Magnetic surfaces 

From the reflective correction for atoms above a di- 
electric surface, the reflective correction for atoms above 
magnetically active surfaces can be found using the du- 
ality of the electric and magnetic fields by applying the 
transformation e <->• /x and d — > — m/c or rh — > cd. For 
example, the near-field reflective correction to the dipole 
spin flip rate for an atom above a magnetic surface can 
be found from Eq. (|5D|) as 



Ik) 



16tt1iz 3 |/i(w) 



{m 2 x + m 2 y + 2m 2 z } . (110) 
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Similarly, the near-field reflective correction to the dipole 
spontaneous decay rate can be found from Eq. (1^7)) as 



r (S) 



H" + dl + 2d 2 z }. (Ill) 



647rfeoc 2 2 

All other results convert in an identical fashion 



other aspect of the far field result is the sinusoidal na- 
ture of the reflective correction with the atom-surface 
distance. This is caused by interference of virtual elec- 
tromagnetic waves reflected from the surface. 



VIII. ACKNOWLEDGEMENTS 



VII. SUMMARY 

Using the theory of electromagnetic field quantization 
in absorbing magnetoelectric materials, the spontaneous 
decay and spin flip rates for the first few multipole or- 
ders have been calculated. In the calculations the primi- 
tive moments, as opposed to the traceless moments, have 
been used since the primitive moments are applicable to 
all physical situations whereas the traceless moments are 
only applicable where the external field obeys the Laplace 
equation and hence is divergence free. The decay rates 
for an atom in free space are given by Eqs. (1551) . (ISUl) . 
(|62|) . (|65|) and §7$, respectively. 

When considering an atom in the prescence of an ab- 
sorbing surface, the free-space decay rate is modified by 
a reflective correction which is a result of wave reflec- 
tion from the surface. In the near-field limit, because 
of absorption by the surface, there are extra degrees of 
freedom to which the atom can decay. The atom can 
lose energy non-radiatively to the surface and hence the 
final state phase space is larger and hence the decay rate 
is increased. The calculation shows that the reflective 
correction is a complicated function of the multipole ori- 
entation but there is also an inverse dependence on the 
atom-surface distance z. The appearence of this type 
of scaling law is consistant with previous studies [L8l — 
[2lj where similar distance relations were found for the 
electric quadrupole decay rate in a variety of systems. 
The ^-dependencies of the dominant order in the reflec- 
tive correction for the different multipole decay rates are 
given in Table ITO1 



surface 


electric 


magnetic 


atom 
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magnetic 


electric 


magnetic 


dipole 


z-* 


z~ l 


z- 1 


z-* 


quadrupole 


z~ 5 


z~ z 


z~ 3 


z~ 5 


octupole 


z~ 7 






z- 7 



TABLE III: Distance dependence of near-field multipole tran- 
sition rates. 

Note that there is a clear hierarchy, with the power 
law dependence of z reducing by 2 with each successive 
multipole order. It is clear from the structure of the mul- 
tipole expansion that this trend will continue to higher 
order multipole moments. 

In the far field the reflective correction is inversely de- 
pendent on the atom-surface distance to all orders of the 
multipole expansion and hence the far fields created by 
all multipoles follow the same distance scaling law. The 
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Appendix A: s-integral parametrization of the 
polarization and magnetization fields 

A common parametrization for the polarization and 
magnetization fields involves re-writing the polarization 
and magnetization fields as an integral over and auxiliary 
parameter s. The derivations of these parametrizations 
are briefly reviewed here. More information can be found 
in Refs. lEHIH- 

A collection of a point charges located at r a can be 
described by the charge distribution 



P(r) = ^2q a S{r - r Q ). 



(Al) 



In most cases, it is convenient to describe the system of 
charges by a coarse-grained distribution 



PaO) = (^2q<}j S(r-r A ) 



(A2) 



The difference between the actual charge distribution 
(|Aip and its coarse-grained average defines the atomic 
polarization field via the implicit relation 



V -Pa(i-) = -p(r)+p A (r). 



(A3) 



Fourier transforming with respect to the spatial variable 
r gives 



ik-PA(k) = y^g a ( 



l-e 



-ik-(r a -r.4) 



(A4) 



This expression can be re- written in terms of an integral 
over an auxiliary parameter s 

P A (k) = V q a e~ lk rA f 1 ds(r a - r A )e~' lk < i -- r ^ s . 

Performing the inverse Fourier transform gives the 
parametrized expression for the polarization field 

Pyi(r) = y2la ds(i a - r A )5 [r-r A - s(f Q - r A )] ■ 

(A6) 
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The magnetization field is generated by currents within 
the material and hence is related to the rate of change 
of the distribution of charges. Differentiating Eqn. (|A3j) 
with respect to time and applying the continuity equation 
p + V • j = gives 



V- j(r)-j A (r)-P A (r) = 0. 



(A7) 



The quantity j(r) — j^i(r) — P A (t) is known as the magne- 
tization current. Since this represents a divergence free 
quantity it can be written as the curl of another vector 
field, the atomic magnetization M^fr) 



V x Mi(r) =j(r) - j A (r) - P A (r). 



(A8) 



For globally neutral ensembles of charges at rest, j J 4(r) 
vanishes. The remaining two terms in Eqn. (|A7I) arc 
Fourier transformed with respect to r. Substituting in 
the parametrized expression for the polarization field, 
performing the time differentiation and taking into ac- 
count the spin contribution to j (r) gives 

ikx]ft A (k)=5^9 a j a e-*-*- 

a 

Jo 

- V q a [ dss(r a - r A )(-tk • * a ) e -*M"i+«(*«-**)] 
+ ^?kx 7 a Se ik '*. (A9) 

a 

Integrating the second term on the rhs by parts and ap- 
plying the vector identity Ax(BxC) = B(A • C) — 
C(A • B) to the remaining terms gives 

Mi(k)=^9 a / dss (r a -r A )xr a 



x e 



-ik-[r A +s(r a -r A )] 



,ik-f Q 



(A10) 



The final step is to perform the inverse Fourier trans- 
form to give the parametrized expression for the atomic 
magnetization field as 



Mx(r) =^2<la / dss (r a -r A )xr a 



X 6 [r - v A - s(r a - r A )] 

+ ^ JuSaS (r - fa) • 



(All) 



Note here that the physical variables have been used. 
For a Hamiltonian, the canonical momenta have to be 
introduced instead of the physical velocities. 



Appendix B: Green function for a planar dielectric 
interface 



The reflective part of the Green function in terms of 
the partial Fourier transform of its components is derived 
in Ref. [46j], the result of which is quoted here. Given 
that 



d 2 k 



(2nf 



R a p(k\\, z, z'u)e 



',.,u ik ir( r n- r li' 



(Bl) 

with r || = (x, y, 0) and k|| = (k x , k y ,0), which are vectors 
restricted to the (x, y)-plane, and feii = |k|||, the compo- 
nents for the reflective part of the Green function for an 
infinitely extended planar dielectric material, which fills 
the lower half space z < 0, are given by 



Rxx 



3 i/3 2>0 (N + |z'|) 



20 



z>0 



JTM 



x TE V 



Pz>0 ,2 r 



R z 



a i Pz>o(\z\ + \ z '\) 



2(3 z >o 



1.2 



(B2) 
(B3) 

(B4) 



f^xy f^yx 



J0 z>o (\z\ + \z'\) 



2/3z>o 



„TM 



~Pz>0 



kx.kii 



k k 



% x"-y TE ^x^y 
2 ^>0 k 2 + r k 2 



(B5) 



20 



z>0 



R - l _p*P*>o(\z\ + W\) 

^Pz>0 



-Pz>ok x 



Rzx j 



„TM 



-/3z>ok y 



-R 



zy- 



(B6) 



(B7) 



The functions r™ and r TE are the Fresnel coefficients 
for TM and TE waves 

r TM _ £ ( UJ )Pz>0 — Pz<0 ^TE _ Pz>0 ~ flz<0 



s(uj)0 z>o + Z<O 1 



Pz>0 + P. 



2<0 



with 



P*<o = y/q 2 e(u)-kl p z >o = ^q 2 - k\ (B9) 
and q = u>/c. 

Appendix C: Rotational Averaging of Tensors 

In many cases physical systems are isotropic in na- 
ture but not explicitly isotropic in their mathematical 
description. Substantial simplification can be achieved 
by averaging of tensorial quantities over the sphere. The 
description of the following method follows closely the 
arguments in Ref. [47j . 



14 



Contraction of a tensor with a rotation matrix whose 
components are the cosines of an angle will rotate the 
coordinate system of the tensor by that angle 



j a i ■ ■ -ot n rp 



(CI) 



By integrating the rotation matrices over the Euler angles 
and dividing by the spherical area the rotational average 
can be found 



with 



1 



rplSO 



dtp / d<p 



sine dei^Z" 



I TOLl---a„\rp 



(C2) 
(C3) 



This is cumbersome even for tensors of low rank. A sim- 
pler method is to contract the tensor with an 'averaging 
tensor'. Owing to its spherical symmetry, any averaging 
tensor of rank n can be written as a sum of the basic 
linearly independent rank-n isomers which in turn can 
be constructed from the two fundamental isotropic ten- 
sors in three dimensions; the Kronecker delta &y and the 
Levi-Civita symbol e^fc. Consider a Cartesian frame of 
reference in three dimensions V with a set of basic rank- 
n isomers where the initial non-isotropic tensor lies 
and a Cartesian frame of reference U with a set of basic 
rank-n isomers in which the rotationally averaged 
tensor lies. The averaging tensor JW is a rank-2n tensor 
that maps between the two spaces, J : V — > V, such that 



Furthermore, 



= . f i 



J(") = g" ■ M • /" , 



(C4) 



(C5) 



where M is a matrix of coefficients. Since the isomers for 
each rank are known all that is required is to determine 
is to find M. Combining Eqs. JC4} and JC5|) yields 

g n = g n ■ M ■ (/" ■ /" ) , 
(9 n -9 n ) = (g n -g n )-M-(r-r). (C6) 

Note that the coordinates of V and U are both Carte- 
sian (merely rotated with respect to each other), hence 
the isomers for the two spaces have the same form and 
magnitude 



(g n -g n ) = (f n -r) = sw 



(C7) 



where is a matrix formed from all possible contrac- 
tions of each isomer in the set. Thus 

S(") = • M • S< n > M = (s("A 1 . (C8) 

Hence, I^") can be determined purely from the isomers. 

The rank-2 averaging tensor, used in the rotational 
averaging of the dipole decay rate, is constructed from 
the single rank-2 isomer = 5 a /3, hence S^ 2 ^ is a single 
number 

S (2) =/ (2) / (2) = ^^ =3) (C9) 

M = 1/3 and thus 



(CIO) 



As a non-trivial example we calculate the averaging 
tensor for a rank-4 tensor used in the rotational averaging 
of the quadrupole decay rate. The rank-4 isomers are 

/l (4) = <W<W H = <W<^ > /J = fiavdpp- (Cll) 

The elements of S*- 4 ) are 

S^=f^f^, (C12) 

so for example 



S[f = /i (4) /i (4) = <W V<W<W = 9, 

= A (4) /i 4) = S aP S^S ail 8p v = 3. (C13) 



Hence. 



S (4) 



9 3 3 



1 



4 -1 -1 



3 9 3 UM=- -1 4 -1 
3 3 9 / 30 \ -1 -1 4 



(C14) 

so that we can write the rank-4 averaging tensor as 



i 



ijkm 









T ^ 








(C15) 

Finally, using the above method, the rank-6 rotational 
averaging tensor, used in the calculation of the octupole 
spontaneous decay rate, is given below. 
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ijkpqr ~ 210 



( SijSkpSgr ^ 
$ij $kq Vpr 
&ij $kr $pq 
^ik^jp^qr 
^ik^jqSp-r 
^ik^jr^pq 
Sip5jfc5q r 
Si p Sj q Skr 
OipO jr®kq 
&iq$jk$pr 
^iq^jp^kr 
$iq5j r 5fzp 
^ir^jk^pq 
Si r 5jp5fcq 

V $ir$jq$kp / 



2 

-5 
2 
2 
2 

-5 
-5 



/ 16 -5 -5 -5 2 2 -5 2 

-5 16 -5 2 -5 2 2 2 

-5 -5 16 2 2 -5 2 -5 

-5 2 2 16 -5 -5 -5 2 

2 -5 2 -5 16 -5 2 -5 

2 2 -5 -5 -5 16 2 2 

-5 2 2 -5 2 2 16 -5 

2 2 -5 2 -5 2 -5 16 -5 

2-5 2 2 2-5-5-5 16 

2 -5 2 2 -5 2 -5 2 2 

2 2 -5 -5 2 2 2 -5 2 

-5 2 2 2 2 -5 2 2 -5 

2 2 -5 2 2 -5 -5 2 2 

2 -5 2 -5 2 2 2 2 -5 

V -5 2 2 2 -5 2 2 -5 2 



2 

-5 
2 
2 

-5 
2 

-5 
2 
2 



2 
2 

-5 
-5 

2 

2 

2 

-5 
2 



-5 2 
2 2 
2 -5 
2 2 
2 2 

-5 -5 
2 -5 
2 2 

-5 2 



16 -5 -5 -5 
-5 16 -5 2 
-5 -5 16 2 
-5 2 2 16 
2-5 2-5 
2 2-5-5 



2 

-5 
2 

-5 
2 
2 
2 
2 

-5 
2 

-5 
2 

-5 
16 

-5 



- 5 \ 
2 

2 

2 

-5 
2 
2 

-5 
2 
2 
2 

-5 
-5 
-5 
16 / 



SasSpriS-zC 

SarjSpjSsC 
SarjSpsSyt; 

SarjSpcd-rS 
0~a(0~f3yO~Sr) 

(C16) 
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